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1 Introduction 

The hierarchy problem of the standard model (SM) can be solved by supersymmetry 
(SUSY). Its minimal realization is the minimal supersymmetric standard model (MSSM) 
and comprehensive reviews are given in p~H3]. In order for the hierarchy problem to 
be solved naturally the SUSY breaking scale and the masses of the super-partners are 
expected to be roughly 1 TeV at most. If this expectation is true, the on-going LHC 
experiment is able to produce the super-partners directly and discover them. What is 
required on the theoretical side is to make both predictions of the various SUSY signals 
and the SM backgrounds as precise as possible. The QCD NLO corrections are generally 
mandatory for the reliable predictions. 

At the modern colliders like the LHC and the future ILC the automatization of the- 
oretical prediction is essentially required. One reason is that the number of modes are 
too large to be treated manually case by case. The other reason is that the QCD NLO 
corrections involve more ingredients than the LO case and especially the corrections to 
the multi-parton processes are too lengthy and too complex to be treated by hand. The 
automatization into a computer code can be realized only by general and practical pro- 
cedures which are executed in some simple algorithmic steps. Here we briefly review the 
publicly available codes for the fixed order perturbative QCD. The LO (tree level) has 
already been automated well in Refs. [4T-T16]. Some of these references include also the 
automation of the phase space integration, the parton shower, the decay, and so on, i.e. 
they are event generators. The NLO real emission corrections are also obtained by the 
same packages because the real corrections enter at tree level. About the NLO virtual 
1-loop corrections, some essential ingredients are automatized and the automation is now 
approaching perfection [71l8l lTTl[T71 - [T9] . Each of the real and virtual corrections possesses 
soft and collinear singularities and they cancel each other for observables [20L22J. For the 
multi-parton processes like 2 — y 4 partons, it is almost impossible to realize the cancelation 
in a fully analytical way. Furthermore the inclusion of a jet algorithm makes the analyti- 
cal phase space integration hopeless. Therefore a general and practical procedure to treat 
such soft /collinear singularities is desirable and was discovered in Refs. [23H3T] . Among 
these references especially the Catani-Seymoure dipole subtraction [23l[2l] is widely used 
for NLO predictions of the multi-parton processes in the so called the 'NLO Wishlist' for 
LHC [32J. During the last two years the NLO corrections to the 2 — y 4 processes in the 
wishlist have been successfully obtained using the dipole subtraction [33TI4T] . In order 
to apply the procedure to the various multi-parton processes, automation is required due 
to the increasing number of the subtraction terms. The algorithm and the formulas are 
sufficiently simple to implement into a computer code. We have automated the procedure 
in the package, AutoDipole, and made it publicly available J42JH3]. The source code is 
written in Mathematica and to obtain the reduced Born matrix elements, the package 
interfaces with the stand-alone version of the MadGraph jU|6] . There are the other four 
packages to automate the subtraction procedure in the frameworks of SHERPA |4"4"] . TeV- 
Jet H5], MadGraph/Event HUE], and HELAC SB]- The FKS subtraction J2S1I2Z1 is also 

1 



automatized in the framework of MadGraph/Event 

Next let us consider the status of the automatic calculation for the SUSY processes, 
typically in the MSSM. The number of fields and interactions in the MSSM are more 
than double in comparison with the SM case, and the number of production and decay 
modes are also much increased. Then the automatic calculation is really needed. It is 
not simple even to calculate the spectrum of the SUSY masses and couplings by hand 
and automation is realized [50H53] . The LO and the real corrections are well automatized 
as in the SM case [54T - I63"] . Event generators also have been constructed [6l fTTll6H473] . 
About the 1-loop virtual correction, the automation has been proceeded by the packages 
for some restricted processes in the colliders [51], for the general processes [66J, and for 
the decay processes [56 |[58l [731476] . About the general treatment of the soft/collinear 
singularities in the QCD NLO corrections of the SUSY processes, the extension of the 
dipole subtraction is constructed [24]. However, there is no available implementation 
of the SUSY cases so far. The purpose of this article is to provide the first publicly 
available package of the dipole subtraction including SUSY processes in the MSSM by 
extending the AutoDipole package. All dipole terms in the procedure include the reduced 
Born matrix elements and AutoDipole interfaces with MadGraph to obtain the matrix 
elements. Since MadGraph includes the MSSM [72], the SUSY extension of AutoDipole 
should be straightforward and we realize it as SuperAutoDipole package in the present 
article. The article is organized in the following way: In Section 2, the extension of the 
dipole subtraction into the MSSM is briefly reviewed with a simple process. In Section 
3, the usage of SuperAutoDipole is explained and some benchmark numbers are shown. 
Section 4 is devoted to the conclusion. 

2 Dipole subtraction in the MSSM 

Let us start this section with the quick review of the MSSM. The gauge groups are the 
same as in the SM. The gauge bosons (spin s = 1) are uniquely embedded in the vector 
multiplets and the super-partners are the gauginos (s = 1/2). The quarks and leptons 
(s = 1/2) are included in the chiral multiplets and the partners are the squarks and 
sleptons (s = 0), respectively. The higgs boson (s = 0) is included in the chiral multiplet 
and the partner is a higgsino (s = 1/2). The sector of Yukawa interaction in the SM, 
which is the source of the fermion masses, is extended to the SUSY version by introducing 
a superpotential. In order to provide masses to both of the up- and down-type quarks, 
a second higgs doublet is required. To avoid baryon and lepton number violation, a 
symmetry called R-parity is imposed, which forbids the proton decay. A small amount 
of R-parity violation is still allowed in the parameter region where the prediction does 
not contradict with the data of the present experiments, as reviewed in Ref. [77] ITS]. 
Hereafter we think only about the MSSM with exact R-parity. Furthermore, explicit 
soft SUSY breaking terms are introduces, typically adding the mass terms of the super- 
partners. The masses of the super-partners are assumed to be roughly at the TeV scale, 



which is the proper scale to solve the hierarchy problem. We here take a closer look at 



(field, super-partner ) 


spin: 


1 1/2 


SU(3) C ,SU(2) L 


(quark, squark) 


(Ql, Ql) 
(Qr, Qr) 


3, 2 

3,1 


(gluon, gluino) 


(9, 9) 


8,1 



Table 1: The spins and the representations in the gauge group SU(3) C and SU{2) L are 
shown. The left-handed quark q L and the super-partner q L represent S , f/(2) i -doublet as 
Ql = (u, (1)l and g^ = (u, d)i, respectively. The right-handed quark q R and the partner 
(Jr are SU (2) L -sing\et as qR = Ur/cIr and qR = Ur/cIr. The quarks and squarks have the 
copy of two more generations as (c, s) and (t, b). 



the colored super-partners. The super-partners of the quark and gluon are squark and 
gluino, respectively, and they are summarized in Table [U The masses are roughly at 
the TeV scale as written above. The fields which possess same quantum numbers can 
be mixed. The mixings between the squarks q~L and qR are large in the case of stop and 
sbottom. The mass eigenstates are written as £1(2) and &i(2), respectively. The mixing 
in the first and second generation can be omitted. The gauge interactions of the squark 
and gluino are derived from the corresponding covariant derivatives, |-D^g| 2 and gp g, 
respectively. The interactions induce the gluon emission from the squark and gluino legs. 
The gluon emission causes the soft/collinear singularity in the real emission correction. 
The appearance of the new source of the singularity in the MSSM requires new dipole 
terms in the dipole subtraction. 

The dipole terms for the massive quarks in the SM and the ones for the squarks and 
gluino typically in the MSSM were constructed at the same time [21]. As the squarks 
have mass but no spin, the dipole term is a bit simpler than the massive quark case and 
the color factor is same as the one of the quark. As the gluino has mass and the spin 
of one half, the kinematical part of the dipole is same as the massive quark case and 
the color factor is same as the one of the gluon. Since the spins of the colored super- 
partners are zero or one half but not one, the SUSY dipoles do not have the correlation 
of the emitter helicity between the amplitude and the conjugate. All explicit formulas 
for the SUSY dipoles are given in [21] . The algorithm to create the SUSY dipoles is also 
straightforward extension from the SM case. Furthermore, since the squarks and gluino 
are sufficiently heavier than the top quark in the MSSM, we do not have to take care 
of the massless limit in practice. Then all dipoles for the SUSY cases are of the same 
type as the gluon emission, which is sometimes called the diagonal splitting type. In the 
AutoDipole package we categorize such diagonal splitting type of dipoles as 'dipole 1' [13]. 
Following the category we add two kinds of the SUSY dipoles to the 'dipole 1' as shown 
at the left in Fig. [TJ For example, the dipole term with the gluon (j) emission from the 
emitter squark(i) in the final state is the case of SUSY (1) at the left in Fig. [TJ If the 
dipole has the spectator(fc) in the final state, it is written as Dy^ = D S quark g iuon,k and 



Dipole term 



SUSY (2) 




I term 



SUSY(l) 




SUSY (2) 



5--'''SUSY(l) 



Figure 1: The emitter of the dipole terms (I-terms) for the SUSY processes are shown in the 
left (right). The symbols q and g represent the squark and gluino respectively. 



the explicit form is written as 

with the dipole splitting function, 

V^fc 



Vy, fc Bl(ij,k) 



(2.1) 
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ijik 
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where Sij = 2pi ■ pj and the quantities, z, y, v, v, are defined in 
color linked Born squared (CLBS) and is denned as 

Bl{ij,k) = (Bl pVTfclBl), 



(2.2) 
The Bl(ij,k) is the 

(2.3) 



where the |B1) represents the amplitude of the reduced Born process in the category 
'dipolel', which is made by the removal of one gluon in the final state from the given real 
emission process. The color operator Tjj(Tfc) acts on the reduced Born amplitude |B1) 
and inserts the color factor into the emitter(spectator) line. The Casimir operators in 
the dipole terms always cancel in the 'dipolel' and they are eliminated in Eq. (12 .ip for 
simplicity. 

The 1-loop virtual correction owns soft and collinear singularities which originate from 
the loop integration and they are regulated as the poles 1/e and 1/e 2 in the dimensional 
regularization. The dipole subtraction procedure provides the I-terms to subtract the 
poles from 1-loop squared matrix element. Corresponding to the SUSY dipole terms two 
kinds of I-terms are added to the SM ones as shown at the right in Fig. [TJ For example 
the I-term with the emitter of the squark(i) and the spectator of a massless parton(k), 
h,k = I S quark,k, belongs to SUSY(l) category. The I-term is explicitly written as 
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i.k 
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2vr T(l - e) C F 
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f^ 
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Sik 



Bl(i,k) 



(2.4) 



with Vi,fc(mi) = VSi t k{mi) + V N Si^ijrii) . Here the singular part VS includes the poles as 

•> , a I \ l * i m i ! i 2 m i ^ 

VSikirrii) = — H In — hr — 

1 ln mf ln s ik 1 ^ m 2 ln s ik 

2 Sjfc Q ife 2 Q ik Q ik 

The non-singular part VNS, and the other quantities, Tq(m.i), jq, Kg are defined in 
Refs. [23j- The CLBSs in a dipole term and in the corresponding I-term are the same 
function of the input momenta. The difference is that the CLBS in the dipole term has as 
input the reduced momenta from the phase space of the real correction process and one 
in the I-term has as input the momenta of the LO phase space. The (Super) AutoDipole 
interfaces with MadGraph to obtain CLBS in four dimension. Then the generated I- 
term by the package belongs to the 't Hooft-Veltman scheme (tHV) in the dimensional 
regularization. The 1-loop correction also must be calculated in the same scheme. But 
the direct application of this scheme to SUSY process causes a problem of SUSY breaking, 
because in the scheme the gluon running in the loop is treated as D-dimensional but the 
super-partner, gluino is treated as 4-dimensional. Fortunately the problem is solved by 
introducing a counter term in such a way as to restore SUSY at the level of the physical 
amplitude in Refs. [79l[80]. About the P- and K-terms their emitter is only a massless 
parton in the initial state and the squarks and gluino in the final state can be only their 
spectator. The kinematical part of the P-term is insensitive to the spin of the spectator 
and the formula of the SM case is not changed for the SUSY extension. Then only the 
K-term must be extended in the case that the spectator is a squark or gluino. Again all 
formulas for the dipole, I-, P- and K- terms are given in [24"] . 

In the remainder of this section, to demonstrate the SUSY cases explicitly we take 
one of the simplest example, the stop pair production, e~e + — > tit*, and will show the 
analytical expressions of the squared matrix element of the real correction |M| 2 cal , the 
dipole term D, and the subtracted real cross section a roal . The squared matrix element of 
the real correction is calculated as 

|M(e^(l)e + (2)^t 1 (3)tt(4)( ? (5))| 2 = 25Qn 3 a 2 e Q 2 a s N c C F 2 \ 2 ■ A (2.6) 

s 12 s 35 s 45 

with 

A = Si2S^S^[Si2{ — S45 + 2^14 + S15) + 335(345 — S14 — S15) + £14(345 — 2Si4 — 2Sis)] 
-2m 2 [Si 2 S 3 5S 4 5 + S 12 (sl 5 (-S 4 5 + Su + S15) + S 3 5S 4 5(-S 4 5 + 2s M + S i5 ) + S45S14) 

— s 35 (si4 + S15) — 2535545^14(314 + S15) — s 45 s 14 ] + 2m 312(335 + S45) (2.7) 

which is uniquely expressed by the minimal set of the independent Lorentz scalars Sij 
and the stop mass m. Here we take the five scalars (si2, S14, S15, S35, S45) as the set. The 
other quantities are written as the electromagnetic coupling constant a e , the charge of 
the stop Q t = 2/3, the strong coupling constant a s , and the usual color factors N c = 3 



and Cf = 4/3. One dipole for the real correction is calculated as 
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Bl(35,4) (2.8) 



with the CLBS, 




Bl(35,4) = -C F - |M(e-(l)e + (2) ^ ^(35)^(4)) | 2 


= -C F N c 87r 2 a 2 e Q 2 t ^- 


s\ 2 - Am 2 s 12 


A(s 2 U5 - Am*) 



(s 34 5 + 2m 2 )(s 34 5 - 2s 35 - 2m 2 ) + s 



:sr, 



Su - 



(3345 ~ S35 + 2m 2 )s 12 
2s 345 + 4m 2 



(2.9) 



where the reduced Born amplitude M(e~e + —¥ tit{) has the reduced momenta P35 and p 4 
for the stop and anti-stop. The S345 is defined as S345 = S34 + S45 + S53. One more dipole 
045,3 is obtained by the trivial replacement of (3 O 4) in Eq. (12. 8p . and the replacements 
of (3 <H- 4) and (1 <H- 2) in Eq. (12. 9ft . The collinear singularities of the squared matrix 
element in Eq. (12. 6 p are screened by the masses of the stop and anti-stop, and the soft 
singularity is subtracted by two dipoles as in Eq. (12. 81) . Then the phase space integration 
of the subtracted quantity becomes finite. In order to see the finiteness concretely we 
here perform the phase space integration in an analytical way. For simpler integration we 
introduce the averaged lepton tensor for the initial electron-positron current and obtain 
the subtracted infrared (IR) safe real cross section as 



a(e~e + -> £~it^) rcal = |d$ (3) ( |M( e - e + -> M^)) 2 - D 35 , 4 - D 34 , f 



a — L-f ■ 1 

7T 



(2.10) 



where the measure of the phase space d<&^ includes the other factors as the fluxes and 
the LO cross section in the massless limit is written as do = nalQ 2 N c /(3si2). The mass- 
and energy- dependent part is calculated as 

VI - 4p 2 



I 



2p 2 



l + Ap- Ap z - lQp 6 + 12 p q 



-2(l-5p 2 + p 4 )ln 



1 - yj - 4p 2 
2p 



-2(l-4p 2 ) 3 / 2 ln^t^ 

P 

i-3p 2 -(i-p 2 )vr^4^ 



V(l-P 2 )ln 



2p 3 



(2.11) 



with p = m/y/si2- Then we can confirm the finiteness of the subtracted cross section with 
help of the analytical expression. In this way the soft /collinear singularities for any real 



emission process in the SM and the MSSM are subtracted by the corresponding dipole 
terms and the subtracted cross section becomes finite. It is noted that for the multi-parton 
processes the analytical phase space integration with a jet algorithm is almost impossible 
and the integration can be performed generally by the Monte Carlo integration. The 
AutoDipole package provides the Fortran routines of the integrand, which can be directly 
integrated by the Monte Carlo method. 

3 Usage of SuperAutoDipole 

The extension of the algorithm for the SUSY cases is straightforward and all the formulas 
are given as seen briefly in the previous section. The calculations of all D-, I-, P- and K- 
terms in the MSSM are reduced to the tree level processes in the same model as shown 
in Eq. (12.11) . In order to obtain the square of the reduced matrix element, our package 
interfaces with MadGraph. MadGraph automatically provides codes of the squared matrix 
elements for the general processes at tree level in the SM [4] . Later it is extended to include 
the processes in MSSM [72] and some models beyond the SM [6]. MadGraph generates all 
diagrams and writes down the Fortran routines to call the corresponding wave functions, 
vertices, and so on in the HELAS library [8TH83] . The choice of the available models is 
well structured and determined at the input file. This feature of MadGraph allows the 
straightforward extension of the AutoDipole package to include the SUSY cases without 
changing the basic structure of the original one for the SM. In the present section we 
would explain the usage step by step with the example in the previous section focusing 
on the new commands. The installation of SuperAutoDipole is exactly same as with 
AutoDipole. The usage is also basically same as 

1. Inclusion of driver : « driver_user.m 

2. Setup of parameters : parameter .m 

interactions . dat 
par am_ card . dat 

3. Run of package : GenerateAll [{initial}, {final}] 

4. Run of the generated codes : make runD, runl , runPK, checkIR 

The steps 1. and 4. are exactly same and here we would explain mainly the new points 
at steps 2. and 3. 



Step 2. Setup of parameters 

parameter. m : 

One new parameter is model whose value can be sm or mssm. If the user input is a 
process in the MSSM, it must be set as model=mssm. The default is model=sm for 
the processes in the SM. One more new parameter is massiveb which specifies the 
masssive(l) or massless(O) of the bottom quark. Default is the massive case. 



interactions.dat : 

The file is one of the MadGraph which is located at . /patch/Models/sm ( . /patch/ 
Models/mssm) for the SM(MSSM) case. User can switch on/off the interactions 
in the SM/MSSM, which is the exactly same way with the use of the MadGraph. 
For the example process we would switch off the Z propagation for simplicity. It is 
realized by the commenting out one line as '#e- e- z GZL QED'. 

param_card.dat : 

The file is also in MadGraph which is at the same place with interactions.dat. 
It specifies the coupling constants, the masses, the decay widths, and so on. If the 
parameter massiveb at parameter. m is set at massless(O), the bottom quark mass 
at this file also should be set at zero. The decay width of a participating heavy 
particle sometime breaks the cancelation of soft/collinear singularities. Then in the 
example, the decay width of the stop should be switched off as 'DECAY 1000006 
. 00000000E+00 # stopl decays'. The parameters of the file can be changed after 
the generation of the code in each created directory under ./process. 

Step 3. Run of SuperAutoDiplole 

The package is run by the same command with the AutoDipole for the SM as 

GenerateAll [{e,ebar},{stl,stlbar,g}] , 

which generates the Fortran routines of all D, I, P, and K terms. The input process 
is the example one, e~e + — > t\t\g. All fields in the SM and the MSSM are available 
and the input forms of the SM and MSSM fields are summarized in Table [2] and 
[3] respectively in Appendix [A] The input forms of the publicly available packages 
can be various, depending on the author's convention, the programming language, 
and so on. In order that the packages can be used by many users comfortably 
without confusion, it is quite reasonable to unify the input form and the definitions 
of the usual parameters. In the community of the automatic calculation of the 
MSSM processes, such necessity is well understood and the efforts have been made 
[MH86] . In the same direction the PDG numbering scheme for the SM, MSSM, and 
some popular extensions of the SM are defined in [EI]. Following the accord, our 
package also provides the input form in the scheme by the command, 

GenerateAllNumber [{11,-11} , -[1000006 , -1000006 ,21}] , 

which is exactly equivalent with the above command. The list of all fields in the 
numbering scheme is also shown in Tables [2] and [31 The command generates the 
Fortran code in the new directory under ./process. The name of the new direc- 
tory follows the convention of MadGraph, which is also shown in Tables [2] and [3j 



The evaluation of the generated code at the example process, e~e + — > t\t\g, com- 
pletely agrees with the analytical result in Eqs. (12. 6p and (12. 8p . In order to provide 
benchmarks for the SUSY dipoles, the values are shown on two phase space points 
in Table H] in Appendix [B] The package also provides the automatic check of the 
process on the user's machine by the command, 

. /lib/sh/check_susy . sh . 

It confirms that the evaluation of the generated code on the user's machine agrees 
with the analytical results in Eqs. (12. 6ft and (12. 8 p on five phase space points, in- 
cluding two used in Appendix [B] When it is confirmed, the command returns the 
message, 

|M|~2(Real) : Confirmed 
Dipole : Confirmed. 

Some advanced users may prefer the run of the package from the shell. It is possible 
by the shell script as 

math « EOF 

<<driver_user . m 

GenerateAll [{e, ebar}, {stl, stlbar, g}] 

Quit 

EOF. 

Such a way is also used at the above mentioned shell script, check_susy. sh. 

The AutoDipole package provides all D-, I-, P- and K- terms, in other words it pro- 
vides all integrands of the phase space integration, except for the 1-loop virtual correction. 
The I-terms are provided by the package in the tHV scheme as mentioned in the Sec. EJ 
On the other hand there is a more suitable scheme for the SUSY process, i.e. dimensional 
reduction (RD) where SUSY is preserved, because the gluon and gluino are treated in 
same dimensionality. The transition rule from tHV to DR is known in Refs. [88H90] . The 
choice of a scheme by one input parameter will be automated in the near future. One 
advantage of the package is that the user can copy the created directory into an user's 
directory and perform the Monte Carlo integration on his favorite setup to obtain the 
observables. The feature easily allows also the user to compare the results of his own 
code with the results of the generated code by the package. Here in order to confirm the 
finiteness of the subtracted cross section as well as the quality of the generated code, we 
execute the VEGAS Monte Carlo integration of the example process on a private setup. 
The 5 x 10 8 phase space points times 10 iteration in the VEGAS algorithm are evaluated 
with the double precision accuracy and the result is obtained in the unit of atto barn as 

a(e~e + ->• iJlg) = 92.5719693 ± 0.0001715 [ab] (3.1) 



where we use the default values of the MadGraph as written in the caption of Table. H] 
and the center-of-mass energy is set at yfs\2 = 4[TeV]. The corresponding value of the 
analytical expression in Eq. ( 12.1 Op is 

a = 92.5719027 [ab] (3.2) 

The results in Eqs. ( 13. ip and (13. 2p completely agree with the sufficient precision. The 
agreement gives a guarantee about the quality of the generated codes. 

Finally we would add the second benchmark numbers for the process, g(l)g(2) — » 
ti(3)ti(4) #(5) which is a sub-partonic process of the real emission correction to the 
observable in a hadron collider, pp — > t\t\ + X. For this process 12 dipoles exist, 
035,4,035,1,035,2,045,3,045,1^ D4 5i 2,Di 5i 3,Di 5i 4,Di 5i 2, D 25 ,3,D 2 5,4, and D 25 ,i. The values 
of the squared matrix |M| 2 and the sum of the dipoles YJ D are shown on two phase space 
points in Table H] in Appendix IE1 The command check_susy.sh includes also the check 
of these values on the user's machine. 

4 Conclusion 

We extend the AutoDipole package to include the SUSY dipoles for the processes in 
the MSSM. The extended package, Super AutoDipole, is the first publicly available one of 
automated dipole subtraction for the SUSY processes. All fields in the SM and the MSSM 
are available. The generated code for the stop pair production is checked against the 
analytical expression including the phase space integration. The package makes possible 
the QCD NLO corrections to the SUSY multi-parton processes like tit\ + jets observables 
at the LHC. Thus the package can effectively help the detailed analysis of SUSY at the 
LHC after they will be really discovered there. The SuperAutoDipole package is available 
for download from [91] or from the author upon request. The MadGraph package (stand- 
alone version) which includes the HELAS library may be obtained from |92j . 
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Appendix 

A Input form 



Field 


PDG 


Input 


Output 


d 


1 


d 


d 


u 


2 


u 


u 


s 


3 


s 


s 


c 


4 


c 


c 


b 


5 


b 


b 


t 


6 


t 


t 


e~ 


11 


c 


e-/e+ 


Ve 


12 


nue 


ve 


^ 


13 


muon 


mu-/mu+ 


^ 


14 


numu 


vm 


T 


15 


tau 


ta-/ta+ 


U T 


16 


nutau 


vt 


9 


21 


g 


g 


7 


22 


gamma 


a 


Z° 


23 


Z 


z 


w + 


24 


Wp/Wm 


w+/w- 



h° 



25 



Table 2: The input forms for the SM fields are shown. The column of 'PDG' shows the 
numbers in the PDG numbering scheme [87]. The numbers are used as the input form 
of the command, GenerateAllNumber [] . The number of the anti-particles is obtained 
by adding minus sign to the particle case. For example, '-1' is used for the anti-down 
quark. The 'Input' shows the input forms of the command, GenerateAll [] . The form of 
the anti-particle is obtained by adding 'bar' to the last of the particle case. For example, 
'dbar' is for d. One exception is the case of W ± and the forms are explicitly shown. 
The 'Output' shows the names of the created directly under ./process. The name for 
the anti-particles is obtained by adding 'x' to the last. For example, 'dx' is for d. The 
exceptions are e~,(j,~,r~, and W + , and shown explicitly. 
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Field 


PDG 


Input 


Output 


d L 


1000001 


sdl 


dl 


UL 


1000002 


sul 


ul 


sl 


1000003 


ssl 


sl 


Cl 


1000004 


scl 


cl 


h 


1000005 


sbl 


bl 


h 


1000006 


stl 


tl 


e L 


1000011 


sel 


el-/el+ 


VeL 


1000012 


snuel 


sve 


P>L 


1000013 


smul 


mul-/mul+ 


V^L 


1000014 


snumul 


svm 


T\ 


1000015 


staul 


tal-/tal+ 


VtL 


1000016 


snutaul 


svt 


d R 


2000001 


sdr 


dr 


UR 


2000002 


sur 


ur 


Sr 


2000003 


ssr 


sr 


cr 


2000004 


scr 


cr 


h 


2000005 


sb2 


b2 


h 


2000006 


st2 


t2 


e R 


2000011 


ser 


er-/er+ 


m 


2000013 


smur 


mur-/mur+ 


T-2 


2000015 


stau2 


ta2-/ta2+ 



1000021 



sg 



N ± 


1000022 


sNl 


nl 


N 2 


1000023 


sN2 


n2 


ct 


1000024 


sClp/sClm 


xl+/xl- 


N 3 


1000025 


sN3 


n3 


N 4 


1000035 


sN4 


n4 


ct 


1000037 


sC2p/sC2m 


x2+/x2 


h° 


25 


h 


hi 


H° 


35 


HO 


h2 


A 


36 


A0 


h3 


H+ 


37 


Hp/Hm 


h+/h- 



Table 3: The input form for the MSSM fields are shown. The notation of the MSSM 
fields in the column of 'Field' follows one in [2]. The definitions of the other three columns 
are same with the SM case in Table. |2J The exceptions are shown explicitly. 
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B Benchmark numbers 

The evaluations of the squared matrix elements and the dipole terms are shown in Table 
H] at two processes, e~e + — > tit\g and gg —¥ t\t\g on two phase space points, PSP1 and 
PSP2 in the unit of [GeV]. 

PSP1 : pi = (500,0,0,500), 

P2 = (500,0,0,-500), 

£>3 = (465.4575514300951,88.15610495656038,197.5104870839467,-100.6674552035491) 

P4 = (442.2757454441043,-9.535905415651449,-180.8891965873525,55.32717033368853) 

p 5 = (92.26670312580057,-78.62019954090894,-16.62129049659421,45.34028486986055) (B.l) 

PSP2 : pi = (500,0,0,500), 

p 2 = (500,0,0,-500), 

p 3 = (443.0578523569805,72.10954684770206,123.4569300473834,-126.9801398447153) 

£>4 = (401.6307218084469,-0.27.27253369182125,-16.78202732992333,23.38599726779297) 

£>5 = (155.3114258345725,-44.83701315588080,-106.6749027174601,103.5941425769223) (B.2) 





|M| 2 [GeV^] 


£D[GeV- 2 ] 


e e + -> tit\g 


PSP 1 


AutoDipole 

Analytics 


3.38015568897016 • KT 7 
3.38015568897014- 10- 7 


4.38408930011336 -10- 7 
4.38408930011337- 10- 7 


PSP 2 


AutoDipole 
Analytics 


7.46080394542307- lO" 8 
7.46080394542308 • 10~ 8 


-2.35778790281811 • 10~ 7 
-2.35778790281812 • 10" 7 


gg -> htlg 


PSP 1 
PSP 2 


2.83987176058353 -10- 4 
2.10737068989263 • 10" 4 


2.77981806065659 • 10~ 4 
2.01941234529644 • 10~ 4 



Table 4: The results of the |M(e~e + — > t\t\g)\ 2 and the sum of two dipole terms J2 D on 
PSP1 and PSP2 are shown. The first row is the evaluation of the automatically generated 
code by the package. The second row is the evaluation of the analytical expressions in 
Eqs. (12.61) and (J27HJ). The results show the perfect agreement of 14 digits. The values of 
the process gg — > t\t\g are also shown at the last two entries. We use the default values 
of the MadGraph as a s = 0.118, a e = 0.00781653039848672 and m = 399.668493[GeV]. 
These numbers can be reproduced on the user's machine by the command check_susy . sh. 
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